NMR study of the dynamics of 3 He impurities in the proposed supersolid state of solid 4 He 
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The dynamics of 3 He atoms in solid 4 He have been investigated by measuring the NMR relaxation times 
Ti , T2 in the region where a significant non-classical rotational inertia fraction (NCRIF) has been reported. For 
3 He concentrations X3 = 16 ppm and 24 ppm, changes are observed for both the spin-lattice relaxation time 
Ti and the spin-spin relaxation time T2 at the temperatures corresponding to the onset of NCRIF and, at lower 
temperatures, to the 3 He- 4 He phase separation. The magnitudes of T\ and Tb at temperatures above the phase 
separation agree roughly with existing theory based on the tunneling of 3 He impurities in the elastic strain field 
due to isotopic mismatch. However, a distinct peak in T\ and a less well-resolved feature in T2 are observed near 
the reported NCRIF onset temperature, in contrast to the temperature-independent relaxation times predicted by 
the tunneling theory. 

PACS numbers: 67.80.bd, 67.30.hm, 76.60.-k 



5h 



1 

C 

O 

o 



> 

On 
O 

m 



X 



The discovery of a non-classical rotational inertia fraction 
(NCRIF) in solid 4 He by Kim and Chan [1, 2] has generated 
enormous interest because the NCRIF could be the signature 
of a supersolid state [3]. Several independent experiments 
[4, 5] have shown that the NCRIF magnitude and tempera- 
ture dependence are strongly dependent on defects such as 
3 He impurities and the quality of the crystals. Recently stud- 
ies of the elastic properties of solid 4 He by Beamish and col- 
leagues [6] have revealed a significant frequency dependent 
change in the elastic shear modulus with an enhanced dissi- 
pation peak having a temperature dependence comparable to 
that observed for the NCRIF. These results suggest that the 
dynamics of the 4 He lattice plays an important role in the low 
temperature bulk properties of solid 4 He and rather than ob- 
serving a phase transition to a supersolid state one may be ob- 
serving a thermally excited dynamical response. It is therefore 
important to study the microscopic dynamics of 3 He impuri- 
ties to better understand their role in the NCRIF and shear- 
modulus phenomena. To meet this need we have measured 
the NMR relaxation times of dilute 3 He impurities in solid 
4 He at low temperature. The NMR relaxation rates are de- 
termined by quantum tunneling (via 3 He- 4 He atom exchange) 
and scattering of the diffusing atoms by the crystal deforma- 
tion field around the 3 He impurities and other lattice defects. 
The NMR relaxation rates are therefore very sensitive to the 
elastic properties of the solid 4 He and on any changes in the 
crystal ground state that would modify the tunneling rate. Al- 
though low-temperature NMR data for 3 He impurities in solid 
4 He have been reported in Refs. 7 and 8 for higher 3 He con- 
centrations X3 > 100 ppm, we report the first NMR data on 
isolated 3 He impurities in the region of x% and T in which 
NCRIF has clearly been observed. (Recently Toda et al. [9] 
reported simultaneous NCRIF and NMR data for a sample 
with £3 = 10 ppm, but the NMR signal was only observable 
in their study from 3 He atoms in phase-separated clusters.) 

The samples were prepared by mixing high purity gases 
and condensing the mixture at high pressure (46.2 bar) into 
a polycarbonate cell that contained a pressure gauge (Fig. 1), 



and then solidifying the samples using the blocked capillary 
method. Thermal contact to the sample was provided by a 
solid silver cold finger extending from a dilution refrigera- 
tor. Details of the cell have been reported elsewhere [8]. A 
first sample with £3 = 16 ppm was annealed for 24 hours 
just below the melting point, while a second sample with 
X3 = 24 ppm was annealed for only 30 minutes. For both 
samples the final pressure measured in situ at low temper- 
ature was 27.75(5) bar, corresponding to a molar volume 
v m ps 20.8 cm 3 . 

Standard pulsed NMR techniques were used to measure the 
nuclear spin relaxation times: magnetization recovery follow- 
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FIG. 1. (Color online) Schematic representation of the low temper- 
ature NMR cell. The preamplifier and tuning capacitor are located 
on a 4 K cold plate located a distance of 1.8 m from the sample cell. 
The RF transmitting and receiving coils are simplified in this figure. 



2 




Temperature (mK) Temperature (mK) 



FIG. 2. (Color online) Temperature dependence of the nuclear spin- 
lattice relaxation time T\ for samples with xz = 16 ppm and 24 ppm 
compared to a high concentration sample (xz = 500 ppm, from Ref. 
8). A peak in Ti is observed for both samples at T ~ 175 mK 
In addition, for the Xz = 16 ppm sample Ti drops by a factor of 
~ 100 below 85 mK due to 3 He- 4 He phase separation. Solid (open) 
circles show data taken on warming (cooling). The dashed green 
lines represent Ti values calculated for the impuriton model [10]. 

ing a spin echo to measure Ti, and a Carr-Purcell-Meiboom- 
Gill [11] multiple-echo sequence to measure T2. Separate or- 
thogonal coils were used for excitation and signal reception to 
minimize sample heating and pulse pickup by the cryogenic 
preamplifier. The coils were isolated from the sample cell 
and allowed to float independently from the sample temper- 
ature which was monitored using a carbon resistor calibrated 
against a 3 He melting curve thermometer. The NMR signals at 
these concentrations are very weak and to improve the signal- 
to-noise ratio we used a low temperature RF preamplifier that 
could operate in a strong applied magnetic field [12]. Even 
with this preamplifier signal averaging up to five days per data 
point was needed to realize an adequate signal-to-noise ratio. 

The studies were carried out for a Larmor frequency of 
wl/27t = 2 MHz as the expected relaxation times extrapo- 
lated from previous studies [13-16] would be prohibitively 
long at higher frequencies. In order to determine the 3 He con- 
centrations of the samples we measured the amplitude of the 
NMR echoes at high temperatures (150 < T < 350 mK) 
which follows the Curie law and compared the signal to a stan- 
dard reference sample with x% = 1000 ppm. 

The observed temperature dependences of the nuclear spin- 
lattice relaxation time T\ for samples with X3 = 16 ppm and 
24 ppm are compared with the dependence for a much higher 
concentration sample (X3 = 500 ppm from Ref. 8) in Fig. 2. 
A pronounced peak at T « 175 mK is observed for both the 
16 ppm and 24 ppm samples in contrast with the weak tem- 



FIG. 3. (Color online) Temperature dependence of the nuclear spin- 
spin relaxation time T2 for samples with xz = 16 ppm and 24 ppm. 
In addition to a strong, hysteretic drop in T2 for xz = 16 ppm below 
100 mK due to phase separation, both samples show poorly-resolved 
features near the temperature at which the Ti peaks are observed, 
T « 175 mK. As T2 reflects the spectral density of 3 He motion near 
zero frequency it is sensitive to slow and static 3 He redistributions 
that do not affect Ti. This may also contribute to the temperature 
dependence of the T2 values observed at T > 200 mK for the xz = 
24 ppm sample, which was cooled after less annealing than the xz = 
16 ppm sample. 



perature dependence observed for the X3 = 500 ppm sam- 
ple. Samples with higher concentrations (500 ppm < X3 < 
2000 ppm) were also studied [8, 17] and had temperature de- 
pendences similar to that of the 500 ppm sample except for the 
expected shift in the phase separation temperature. The posi- 
tion of the peak in T\ corresponds closely to the saturation 
temperature T90 observed for NCRIF at this x$ [4]. Allow- 
ing for expected shifts with pressure and X3, the T\ peak we 
observe also correlates well with the shear dissipation peak 
observed by Syschenko et al. [6]. Finally, the T\ peak occurs 
at approximately the same temperature for this X3 value as the 
sharp ultrasonic absorption anomaly reported by by Ho et al. 
[18]. 

In parallel with the observations for Ti, a less well resolved 
peak in Ti is observed for X3 = 16 ppm at the same tempera- 
ture for which the T\ peak is seen (Fig. 3). At lower tempera- 
tures large changes in Ti and T 2 are observed due to 3 He- 4 He 
phase separation at T ps w 85 mK. This value of T ps agrees 
well with the predictions of Edwards and Balibar [19]. The 
same features are observed for the sample with 24 ppm: a 
strong peak in T\ at T « 175 mK (Fig. 2) and a weaker peak 
in T 2 (Fig. 3). The Ti and T 2 features at T w 175 mK ap- 
pear unrelated to the large changes in Ti and T 2 associated 
with the phase separation at lower temperatures. To verify 
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this fact, the 24 ppm sample was purposely never cooled to 
the phase-separation temperature. 

The thermal hysteresis around 85 mK in T\ and T2 rein- 
forces the view that the changes at this temperature are due 
to 3 He- 4 He phase separation [20]. The NMR echo ampli- 
tudes become temperature independent below the phase sepa- 
ration as expected for the formation of degenerate liquid 3 He 
droplets [17]. Our results for the phase separation tempera- 
tures are in good agreement with those reported by Toda et al. 
forxa > 100 ppm [7]. 

The samples with x% = 16 ppm and x% = 24 ppm have 
significantly different temperature dependences for T2 at T > 
200 mK. As T2 (unlike Ti) is sensitive to very slow and static 
changes in the positions of the 3 He impurities, it seems that 
the T2 data cannot clearly resolve the effects leading to the T\ 
peak from other effects. In particular, the higher-temperature 
data in Fig. 3 suggest that T2 is sensitive to crystal quality. 

The nuclear spin dynamics of 3 He impurities in solid 4 He 
have been studied extensively for relatively high concentra- 
tions (X3 > 90 ppm) and high temperatures (T > 350 mK) 
[13-16]. The results have been described in terms of mobile 
3 He impurities tunneling through the 4 He matrix by 3 He- 4 He 
exchange ( J34) with a mutual scattering due to the elastic de- 
formation field surrounding each impurity (K(r) = J^r -3 ) 
[10, 21]. A reasonable fit to the I\ and T 2 values observed 
for very low concentrations (~ 20 ppm) at high temperature 
(T > 200 mK) is obtained using the Landesman model [10] 
as shown by dashed green lines in Figs. 2 and 3. The 3 He 
atoms can also tunnel as weakly bound pairs [22]. The pair- 
tunneling model explains resonant dips in T\ observed [15] 
at = 1-3 and 2.6 MHz since T\ directly measures fluctu- 
ations in 3 He- 4 He separations at frequencies of wl and 2ojl- 
However, none of these models predict temperature depen- 
dent relaxation times so they can not explain the observed T\ 
peaks. 

The T\ peak occurs at roughly the same temperature at 
which torsional oscillator and shear modulus anomalies are 
observed [4, 6], and it is tempting to infer a connection be- 
tween all three phenomena. One possibility is that the sharp 
reentrant peak in fluctuations near 175 mK signals a phase 
transition, possibly associated with supersolidity. However, 
other explanations not involving a phase transition must be 
considered. As shown in Ref. 4, the £3 -dependent onset tem- 
perature for the torsional-oscillator anomalies agree well with 
the temperature Tip(x3) below which 3 He pinning of disloca- 
tion lines is expected to dominate over pinning by dislocation- 
network nodes, Tip = B / \x\(a / x^L-^) where B w 0.5 K is 
the binding energy of a 3 He impurity to a dislocation line, 
i N is the average length of dislocation segments between 
nodes, and a = 3.7 x 10~ 10 m is the nearest-neighbor dis- 
tance. In this model there are roughly L N /a binding sites 
for a 3 He impurity on an internode segment of a disloca- 
tion line, and each binding site is occupied with probability 
P oc = (e- B / T /x 3 + l)-\ 

For a total length of dislocation lines per unit volume A 
in 



capillary method [4], the concentration of 3 He binding sites on 
dislocations relative to the number of 4 He atoms in the sample 
is xa ~ 10~ 8 , much smaller than the concentrations of 3 He 
atoms X3 present in our NMR experiments. It would be sur- 
prising if such a small concentration of binding sites Xd could 
have a measurable effect on T\ . However, there are other in- 
dications that the density of dislocations or other defects must 
be much larger than the estimate quoted above, if the NCRIF 
observations are explained either by superflow along disloca- 
tions cores [23] or as a direct mechanical effect unconnected 
with supersolidity [24]. 

Figure 4 shows the temperature dependence of T\ along 
with the occupation probability for 3 He binding sites and 
another model discussed below. With the expected 3 He- 
dislocation binding energy B « 0.5 K it can be seen that the 
significant transition in P oc (T) occurs at much lower temper- 
atures than the peak in T\ as shown in Fig. 4(b). If larger 
binding energies B = 1.8 — 1.9 K are assumed, the drop 
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FIG. 4. (Color online) (a) Measured spin-lattice relaxation times for 
two low-a;3 samples (data points) along with fits to a phenomenolog- 
ical model based on thermally-activated relaxation of unknown de- 
grees of freedom (smooth curves), (b) Temperature-dependent prob- 
ability for a defect such as a dislocation to be occupied by a 3 He 
impurity. Here B is the 3 He thermal activation energy (a) or defect 



0.2/ w 10 11 m~ 2 for 4 He crystals grown by the blocked- binding energy (b). 
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in P 0C (T) occurs at temperatures close to those of the NMR 
anomaly. However, it is not clear how the partial occupation 
of 3 He binding sites would lead to reduced fluctuations of 3 He 
interatomic vectors as inferred from the peak in T\. 

Another phenomenological model we consider is to asso- 
ciate the T\ anomalies with a thermally activated relaxation 
peak as has been used successfully to describe the shear- 
modulus shifts [6]. The smooth curves in Fig. 4(a) are fits to 
theformTi = [R -2R!Ujt/(1 + (ujt) 2 )}-\ lot = ujT e B / T . 
Here Rq and R\ are fitting parameters giving the background 
relaxation rate and height of the peak in T\ , B is the activation 
energy, tq is an attempt frequency, and u> is the frequency at 
which relaxation is being probed (e.g. uj l = 1.26 x 10 7 s _1 ). 
For the dotted and dashed curves the activation energy was 
fixed at B = 0.5 K (the approximate binding energy of a 
3 He impurity to a dislocation inferred in Ref. [4]) while for 
the other curves B was allowed to increase to improve the 
fit to the data. Note that these fits determine only the com- 
bination lutq (not u) or to separately); for the best fit to the 
x 3 = 16 ppm data cjt = 2.9 x 1(T 5 and B = 1.81 K. 

It must be emphasized that the fitting functions used in Fig. 
4(a) are phenomenological and not derived from a simple mi- 
croscopic model of the NMR relaxation. For example, a sim- 
ple resonance between the thermally-activated motion of 3 He 
impurities and the Larmor frequency wl would ordinarily lead 
to a minimum in T\ (increased fluctuations at the Larmor fre- 
quency), not a peak as we observe. On the other hand, relax- 
ation of 3 He impurities in solid 4 He is known to be controlled 
by a nonmonotonic spectral density of fluctuations with a 
sharp feature at the frequency (~ 3 MHz) at which nearest- 
neighbor 3 He pairs "walk" through the lattice via quantum 
tunneling [22]. Therefore a more complex mechanism such 
as the disruption of the quantum walking motion of 3 He pairs 
by resonant fluctuations of dislocation lines might be needed 
to explain our NMR data. 

The NMR measurements can only be understood in terms 
of a sharp change in the fluctuation spectrum and the same 
changes would also play a dominant role in the anomalies 
observed for the sound attenuation and the shear modulus. 
Specifically, if the fluctuation spectrum is associated with crit- 
ical behavior at a phase transition, the attenuation and NMR 
relaxation rates would vary as T\ ~ |T — T | A in agree- 
ment with the NMR and sound attenuation results. Alterna- 
tively, if a collective but non-critical change occurred in the 
lattice dynamics, the associated change in the elastic proper- 
ties of the solid would result in changes of both the observed 
shear modulus and the NMR relaxation rates. The latter has 
been shown to depend on the elastic strain surrounding 3 He 
impurities [10]. Further studies for a wide frequency range 
and for samples grown at constant pressure to produce higher 
quality crystals are needed to distinguish between these inter- 
pretations. However, it is clear that a successful microscopic 



model of the lattice dynamics of solid 4 He must explain not 
only the NCRIF observations, but also the coincident anoma- 
lies that have been observed in ultrasound [18], shear modulus 
[6], and NMR relaxation as reported here. 
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